1  gt;  BggQ3BHBHBlB32E 


LI=U«I;IW»T 


uvmm 


1.  Rf  port  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


3a  S£Cu«'Tv  CLASSIFICATION  AUTHORITY 


3b  OCCLASSIFICATION/OOWNGAAOING  SCHEDULE 


«  PERFORMING  ORGANIZATION  REPORT  NUMSERISI 


te  NAME  OF  PERFORMING  ORGANIZATION 


■puiiM4.isiu.Mi 


tlf  epplictble  i 


Ecodynamics  Research  Associates 


Sc.  ADDRESS  ICity.  Slat*  and  ZIP  Coda) 

Albuququerque ,  NM  87198 


lb.  RESTRICTIVE  MARKINGS 


3  OISTRIBUTION/AVAILABILITV  of  report 

Approved  for  public  release;  distribution 
unlimited. 


S.  MONITORING  ORGANIZATION  REPORT  NUMSERISI 

AFOSR-TR-  3  3-  0  708 


7a.  NAME  OF  MONITORING  ORGANIZATION 

Air  Force  Office  of  Scientific  Research 


7b.  ADORESS  ( City.  Slab  and  ZIP  Cod*! 

Directorate  of  Mathematical  &  Information 
Sciences,  Bolling  AFB  DC  20332-6448 


Ba.  NAME  OF  FUNDING/SPONSORING 
ORGANIZATION 

AFOSR 


Sc.  AOORESS  ICity.  State  and  ZIP  Code) 


ORCflSE  SYMBOL  B.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
Of  eopUceble) 

9M  F49620-84-C-0079 


10.  SOURCE  OF  FUNDING  NOS. 


PROGRAM 
ELEMENT  NO. 

61102F 


PROJECT 

TASK 

NO. 

NO. 

2304 

A3 

WORK  UNIT 
NO. 


Bolling  AFB  DC  20332-6448 


11.  TITLE  Unclad*  Security  CUmiftcotion) 

Adaptive  Grid  Generation  Using  Elliptic  Generating  Equations  With  Precise  Coordinate  Contro 


13.  PERSONAL  AUTHORIS) 


1 13b.  TIME  COVE 

RED 

1 

14.  DATE  OF  REPORT  |Yr..  Mo..  Dmy) 

wmsimm 

HE 

5  31  July  1985 

COSAT  I  CODES 


GROUP  I  SUB.  GR. 


IB.  SUBJECT  TERMS  (Continue  on  reverte  if  ntctavy  and  identify  by  Block  number) 

symbol  manipulation  codes,  variational  grid  generation 


IB.  ABSTRACT  (Continue  on  reverte  if  neceteery  and  identify  By  Block  number/  <0. —  ~ 

The  investigators  are  in  the  final  stages  of  developing  a  "toolkit*  of  symbol  manipulation 
codes  for  variational  grid  generation.  They  will  present  an  invited  paper  on  this 
work  at  an  AIAA  Aerospace  Sciences  Meeting.  They  havef-ialso  discovered  an  unexpected 
folding  of  the  grid  for  a  design  case  using  a  popular  grid  generation  code.  A  paper 
analyzing  the  problem  was  presented  at  the  AIAA  Computational  Fluid  Dynamics  Meeting 
in  July.  The  analysis  has  suggested  practical  values  of  relative  weighting  parameters 
for  use  in  variational  grid  generation  techniques.  A  paper  has  been  submitted  to  the 
journal  Applied  Mathematics  and  Computation.  ' 

)'  -'-'-t&Zk  DTin 

_  _ •»  rvAnV  C  a  /-A  ELEC-. 


JfflC  FILE  COPY 


SEP  13 1985 


20.  OlBTRIBUTION/AVAILABILITY  of  abstract 
UNCLABBIPISO/UNLIMITEO  E  SAMS  AS  RPT.  □  OTIC  USERS  □ 


no.  name  op  responsible  inoivioual 

3ohn  P.  Thomas 


OD  FORM  1473,83 


31.  abstract  security  CLASSIFICA 
UNCLASSIFIED 


33b.  TELEPHONE  NUMBER  1 23c  OFFICE  SYMBOL 


22b.  TELEPHONE  NUMBER 
(fnctude  Atom  Code/ 

(202)  767-  5026 


unr.iJi.i-MM.'ii.ii'T?'  i 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE 


A 


s 


ECODYNAMICS 
Research  Assoc i at es , Inc . 

P.  0.  Box  B172 
A1 buquqerque ,  N.M.  B7198 
U.S. A. 

(505)  262  -  0440 
31  July  1985 


Captain  John  F.  Thomas.  Jr. 

AFOSR 

Directorate  of  Math,  &  In-fo.  Sciences 

Building  410 

Boiling  AFB,  DC  20332 

Dr.  James  D.  Wilson 
AFOSR /NA 

Directorate  o-f  Math.  &  Into.  Sciences 

Building  410 

Bolling  AFB,  DC  20332 


COPY 

•NSPLCTE0 


Accession  For 

NTIS  GRA&I  ’ 
DTIC  TAB 
Unannornc^-d 
Justification. 


By _ _ 

_Distrlbutl#n/ 

_ Availability  Codes 

j Avail  and/or 
Dlst  j  Special 


SUBJECT:  Research  Progress  and  Forecast  Report  tor  Contract 
F49620-84-C-0079 ,  "Adaptive  Snd  Generation  Using  Elliptic 
Generating  Equations  with  Precise  Coordinate  Controls" 

Dear  Captain  Thomas  and  Dr.  Wilson: 

This  letter  constitutes  the  required  report  on  the  sub¬ 
ject  contract.  Generally,  progress  continues  to  be  excellent. 
The  topical  headings  below  -follow  those  in  our  proposal. 

1.  Use  o-f  Symbolic  Mani  pul  at  i  on . 

The  use  o-f  computer  Symbolic  Manipulation  to  perform  the 
theoretics]  manipulations  for  the  several  variational  -formula¬ 
tions  and  to  produce  the  Fortran  subroutines  tor  numerical  wort 
continues  to  be  productive.  We  are  in  the  Final  stages  c+ 
developing  a  "toolkit"  o-f  Symbol  Manipulation  codes  For  varia¬ 
tional  grid  generation,  and  we  will  be  presenting  an  invites 
paper  on  this  work  at  an  ArtiFicial  Intelligence  session  at  the 
AIAA  Aerospace  Sciences  Meeting  this  January  in  Reno.  (This 
work  is  also  partially  sponsored  by  AR0. ) 

2.  Variational  Formulation  and  ReFerence  Grid  Definition. 

A  summary  paper  on  variational  grid  generation  prin¬ 
ciples,  with  several  variational  formulations,  has  been  prepared 
and  submitted  For  publication. 


Variational  Formulation 
Physical  Space. 


Logical  Space  rather  than 


Approved  for  public  release 
distribution  unlimited. 


This  -fundamental  departure  -from  the  historically  popular 
approach.  as  described  previously.  We  have  now  re-interpreted 
the  "smoothness"  operator  in  logical  space,  and  -find  that  it 
really  is  more  meaning-ful  to  view  it  as  a  constant  segment 
length  condition.  This  i nterpretat i on  clearly  indicates  whv 
-folding  occurs  with  this  -formulation  and  the  more  common 
smoothness  in  physical  space,  and  further  suggests  a  cure  +or 
the  folding  which  involves  the  judicious  use  of  our  “reference 
grid"  approach.  The  key  is  that  the  segment  length  speci fiction 
should  vary  with  the  region  boundar 1 es , 1 n  an  easily  prescribec 
way.  We  are  expecting  this  approach  to  significantly  improve 
the  grids  obtained  with  logical  space  "smootting"  variational 
principles;  testing  and  validating  of  this  concept  and  t~e  ether 
varitional  formulations  will  occur  in  the  next  reporting  period. 

Our  work  on  the  unexpected  folding  of  the  Winslow  (homo¬ 
geneous  Thompson-Thames-Mast 1 n )  method  (which  supposedly  alwa/s 
gives  a  non-fclded  grid  generation  in  2D  because  of  the  max imum 
principle  of  the  homogeneous  equations)  was  receives  w. tn 
considerable  interest  at  the  AIAA  Computational  Fluid  Dynamics 
Meeting  in  July.  The  analysis  was  based  on  the  closed  fern 
solutions  for  minimal  grid  resolutions,  i.e.  a  3x3  gric.  V e 
have  applied  this  approach  more  extensively  now,  to  analyze  nc: 
on l y  the  folding  of  the  Winslow  method,  but  also  the  folding  cf 
•/ar  i  at  i  or  al  1  y  produced  grid  generation  equations.  The  analysis 
has  suggested  practical  values  of  relative  weighting  parameters 
for  the  various  variational  eouations,  and  numerical  experimer- 
cation  in  finer  grids  indicates  that  the  values  so  obtained  are 
indsec  practical.  A  summary  paper  on  these  results  has  been 
submitted  far  publication.  This  work  has  also  shown  that 
variational  volume  control  alone  does  in  fact  produce  a  unique 
system  of  eouations  in  the  minimal  and  resolution,  contrary  to 
the  theorem  cited  b>  Brackbi  1  1  and  Saltzman.  We  are  worur.g  or 
a  correct  interpretation  of  these  conflicting  results. 

4.  Solution  Procedure  Minimizing  Storage  Requirements 

and/or  Arithmetic  Operations 

Tt  i  s  st-atecy  involves  a  solution  formulation  ir  wn.cr 
the  matrix  equation  is  expressed  as  a  variables  sepa'  adie 
elliptic  operator,  with  a  nonhomogeneous  term  which  includes  all 
the  non  1 1 near i 1 1 es  and  non-separat 1 e  coefficients.  Th o 
performance  has  been  excellent,  provided  that  the  grid  genera¬ 
tion  parameters  are  such.  as  to  generate  a  good  grid.  Wt  have 
incorporated  this  formulation  into  a  single  equation  multi  grid 
(FAS)  algorithm,  with  improved  convergence  rates  tor  tine  x.esr 
cal cui at i ons . 

This  work  led  to  a  consideration  of  sub-gric  and 
super-grid  coefficient  generation  for  "black  box"  multigrid 
codes  which  has  been  fruitful.  A  concept  evolved  for  generation 
of  sub-arid  coefficients  given  only  the  discretization  cn  the 
fine  arid,  without  reference  to  or  knowledge  of  the  underlying 
continuum  equations.  The  concept  was  further  extended  in  the 
opposite  direction,  allowing  primitive  coefficient  generatio-  on 
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something  less  than  the  -finest  grid  o-f  interest,  followed  by  an 
economical  ‘‘prolongation"  to  the  finer  grids.  The  concept  is 
most  applicable  to  systems  wherein  the  cost  per  node  point  of 
generating  the  coefficients  is  high,  precisely  the  case  with  the 
variational  formulations.  It  is  expected  to  result  in  a  savings 
by  a  factor  of  5  (limit  8)  for  3D  grid  generation.  The  work  was 
presented  at  the  Second  Copper  Mountain  Multi  grid  Conference  in 
March,  and  interest  was  high  from  attendees  including  A.  Brandt , 
J.  Dendy,  and  U.  Trottenberg. 

5.  Collaborative  Work  with  Dr.  Eisemar.. 

Collaboration  with  Dr.  Eisenan  continues  to  be  quite 
helpful.  We  have  enjoyed  discussions  with  him  prior  to  and 
during  the  AIAA  CFD  Meeting  in  Cincinnati  this  July. 
Discussions  have  covered  a  wide  range  of  topics  related  to  grid 
gene-ation;  particularly  noteworthy  are  the  topics  of  the 
uniqueness  of  the  varitional  volume  control  equations  (see  MZ 
above',  surface  and  curvature  variational  formulations,  compu¬ 
ting  and  graphics  needs  and  capabilities,  behavioural  erro-s 
such  as  grid  folding  and  jacobian  positivity,  evaluation  of 

competing  grid  generation  concepts,  judging  grid  quality,  ana 

his  own  approach  of  iterative  re-di st r l but i on  of  segment 
lengths.  We  will  be  making  a  trip  to  visit  him  at  Columbia  this 

spring,  and  will  provide  him  with  our  codes  and  assistance  in 

installing  them  at  his  facility. 

6.  Interest  from  Air  Force  Weapons  Laboratory. 

We  have  not  heard  any  more  from  AFWL  personnel  here  i  r. 
Albuqueraue  who  had  previously  expressed  interest  in  contracting 
to  create  a  PC  version  o'  some  cf  our  grid  generation  methods. 
We  will  be  meeting  in  mid-August  to  discuss  future  funding  for 
our  continued  development  of  the  ELF  electrode  design  codes,  ana 
we  will  request  funds  to  incorporate  our  new  variational  methodr, 
and  multigrid  techniques  (developed  under  the  subject  contract; 
into  the  2D  and  3D  ELF  codes. 

2.  Publications  Partially  Supported  Under  the  Freser.t 
Contract . 

7.1  "Symbolic  Manipulation  and  Computational  Fluic  Dynamics  , 

Stanlv  Steinberc  and  Patrick  J.  Roache,  Jour.  Compu.t  at  1  or  a  1 _ 

Physics,  Vol  .  57,  No.  2,  Jan.  1985,  pp.  251-284. 

7.2  “  App  1  .  cat  l  or.  of  a  Si  ngl  e-Eouat  i  on  M3-FAS  Solver  to  Elliptic 
Grid  Generation  Equations  (Sub-grid  and  Super-grid  Coefficient 
Generation)’ ,  Patrick  J.  Roache  and  Stanly  Steinberg,  Proc. 
Second  Copper  Mountain  Conference  on  Multi  grid  Methods,  1-3 
April  1985,  Copper  Mountain,  Colorado.  To  appear. 

7.3  "Variational  Grid  Generation",  Stanly  Steinberg  and  Patrick. 

J.  Roache,  submitted  for  publication  to  Numerical  Methocs _ *  or 

Partial  Differential  Equat i ons . 

7.4  "A  Tool  Kit  of  Symbolic  Manipulation  Programs  for  Varia 
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tional  Grid  Generation",  Stanly  Steinberg  and  Patrick  J. 
Roache,  AIAA  Aerospace  Sciences  Meeting,  6-9  Jan.  1986,  Rene, 
Nevada.  To  appear. 

7.5  "A  New  Approach  to  Grid  Generation  Using  a  Variational 
Formulation",  Patrick  J.  Roache  and  Stanly  Steinberg,  AIAA 
-S5-1527-CP,  Proc.  AIAA  7th  Computational  Fluid  Dynamics  Confer¬ 
ence,  15-17  July  1985,  Cincinnati,  Ohio. 

7.6.  "On  the  Folding  of  Numerically  Generated  Grids",  Jose 
Castillo,  Stanly  Steinberg  and  Patrick  J.  Roache,  submitted  to 
app1  le  d  _M  at  hem  at  i  cs  and  Computation . 

7.7  "Electric  Field  Calculations  Using  the  ELF  Codes",  M.  .on 
Dadelszen,  W.M.  Moeny  and  Patrick  J.  Roache,  Proc.  IEEE  F'uisec 
Power  Conference,  Crystal  City,  DC,  10-12  June  1985.  To  appea'. 

7.8  "The  ELF  Codes:  Electrode  Design  for  Lasers  and  Switches', 

F'atrici-:  J.  Roache,  Invited  Paper,  Proc.  CT  AC-85  Conference, 

Melbourne,  Australia,  25-28  August  1985.  To  appear. 

8.  Presentations  Partially  Supported  Unde^  the  Present 
Contract . 

8.1  "Symbolic  Manipulation  and  Computational  Fluid  Dynamics  , 
Patrick  J,  Roache,  Mechanical  Engineering  Seminar,  University  cf 
California  at  Davis,  2  May  1985. 

8.2  "Symbolic  Manipul  at  i  or.  and  Computational  Fluid  Dynam:  cs"  - 
Patrick  J.  Roache,  Applied  Mathematics  Seminar,  Sane:?.  Mat!  oral 
Laboratories  Livermore,  Livermore,  California,  7  Ms v  1935. 

8.3  "Variational  Grid  Generation',  Stanly  Steinberg,  Applies 
Mathematics  Seminar,  Sandia  National  Labo^ator i es  Livermore. 
Livermore,  California,  14  May  1985. 

8.4  "A  New  Approach  to  Grid  Generation  Using  a  Variational 
Formal  at  l  or. "  ,  Patrick  J.  Roache  and  Stanly  Steinberg.  AIAA  th 
Computational  Fluid  Dynamics  Conference,  15-17  July  19E5. 
Cincinnati,  Ohio. 

(future  pr esentati ons) 

8.5  "The  ELF  Codes:  Electrode  Design  +  or  _asers  anc  Switches'', 
kevnote  Lecture,  Patrick  J.  Roache,  Coit.dl  ta  t  -  oral  Techniques  anc 
Applications  Conference,  Royal  Melbourne  Institute  of  Techno¬ 
logy,  Melbourne,  Australia,  25  28  August  1985. 


8.6  "A  Tool  tut  of  Symbolic  Manipulation  Programs  for  Varia¬ 
tional  Grid  Generation",  Invited  Paper,  Starlv  Steinberg,  AIAA 
Aerospace  Sciences  Meeting,  6-9  Jan.  1986,  Reno,  Nevada. 


9. 


Design  Optimization 


In  regard  to  my  -future  work  on  design  optimization 
utilizing  the  ELF  codes,  my  visit  to  AFOSR  on  IS  July  was  most 
helpful.  The  pre-publication  paper  of  Polak's  providec  by  Dr. 
Mark  Jacobs  proved  to  be  just  what  I  need.  I  have  contacted 
Polak  and  he  has  been  receptive  to  investigating  the  use  of  hi s 
algorithms  for  the  electrode  design  problem.  We  had  scheduled  a 
visit  tor  me  to  Berkeley  in  August,  but  I  will  have  to  cancel 
because  of  the  press  of  final  report  preparation  on  another 
contract,  and  my  visit  to  Australia.  I  may  not  be  able  to  visit 
him  until  October  (due  to  jury  duty  in  September)  but  I  defi¬ 
nitely  will  pursue  this. 

10.  Computer  System  Purchase 

Vernita  Slater  of  AFOSR  says  that  I  may  expect  the 
computer  money  in  late  August.  I  very  much  appreciate  your  vote 
of  confidence. 

I  did  not  receive  the  requested  funding  from  ARC  for  the 
rest  of  the  system  I  hac  described  in  my  Request  for  Additional 
Equipment  Purchase  of  7  February  1965.  The  mathematics  people 
at  ARC  were  willing,  but  ARC  policy  prohibited  it  for  legal 
reasons.  ARC*  is  willing  to  charge  the  contract  for  the 
increased  operating  costs  if  Ecodynamics  would  purchase  the  rest 
o-  the  system,  but  this  will  cause  great  accounting 
difficulties,  and  we  would  likely  run  afoul  of  DCAA  auditing 
r egul at : ons . 

Further,  the  ccst  of  the  MicroVa:;  II  is  not  as  attractive 
as  1  har  teen  led  to  believe.  The  hardware  costs  are  close  tc 
the  estimates  1  had  been  given,  but  the  software  costs  are 
out i ardi sn .  DEC  wants  #7K  (after  discount)  for  the  VMS 
operating  system  and  documentation,  and  another  *7F  for  UNIX. 

by  contrast,  the  SUN  microsystem  (which  Prof.  Steinberg 
is  now  getting  up  at  the  university)  is  lower  in  cost  f  o-*  the 
harcware.  and  includes  UNIX  (Berkeley  4.2),  Fortran  77,  Pascal, 
"C‘  ,  and  Core  Graphics  Software  for  no  ewtra  cost.  The 
discounted  price  for  the  SUN  is  only  slightly  over  the  121.5!’ 
requested  from.  AFD3R ,  and  Ecodynamics  can  make  up  the 
difference.  The  present  SUN  based  on  the  virtual  memory  16/32 
bit  MC  68010  is  faster  than  the  MicroVa:;  I,  and  the  upgrade  to 
the  full  32  bit  MC  68020  (available  in  November)  is  eicpectet  tc 
be  somewhat  faster  than  the  MicroVa;-;  II.  Skye  systems  :s 
bringing  out  their  Warrior  array  processor,  a  15  Megaf 1  op  device 
for  $  1 2K ,  including  a  1024  FFT ;  their  two  primary  ta-get 
computers  a^e  the  MicroVa;:  II  and  the  SUN.  (I  would  seel 
additional  funding,  probably  tr0m  AFWL,  for  the  array 
processor.)  I  have  also  investigated  (with  consulting  from 
Professional  Computer  Consultants  of  Albuquerque)  two  alterna¬ 
tive  systems,  tne  ATS<T  UNIX  7300  and  the  IBM  AT  with  the  DSI-32 
co-processor .  Although  both  are  cheaper  than  the  SUN,  the  SUN 
appears  to  be  preferable. 

The  only  advantage  le+t  to  the  VAX  is  compatibility  to 
the  VAX  world,  which  may  not  be  worth  the  price.  Also,  the  SUN 


compatibility  with  Steinberg's  system  would  be  ar  advantage. 
Note  also  that  Prof.  Polak  at  Berkeley  has  a  SUN,  Los  Alamos 
Labs  have  adopted  SUN  as  their  standard  workstation,  ana  Dr.  N. 
L.  Rappagnani  at  AFWL  uses  a  MC  68000  machine  with  Berkeley  4.? 
UNIX  -for  his  grid  generation  work. 

1+  Steinberg's  experience  is  positive,  I  will  probably  te 
requesting  a  change  in  my  previous  Request  For  Additional  Eoui p- 
ment  Purchase  from  the  MicroVax  II  to  the  SUN  System  2/130. 


Thank  you  both  for  your  conti nuec  support. 


